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Abstract 

Using the fluctuation exchange approximation and a three-orbital model, we study the band 
renormalization, Fermi surface reconstruction and the superconducting pairing symmetry in the 
newly-discovered iron-based superconductors. We find that the inter-orbital spin fluctuations lead 
to the strong anisotropic band renormalization and the renormalization is orbital dependent. As a 
result, the topology of Fermi surface displays distinct variation with doping from the electron type 
to the hole type, which is consistent with the recent experiments. This shows that the Coulomb 
interactions will have a strong effect on the band renormalization and the topology of the electron 
Fermi pocket. In addition, the pairing state mediated by the inter-orbital spin fluctuation is of an 
extended s-wave symmetry. 

PACS numbers: 74.70.-b, 74.25.Jb, 71.18.+y, 74.20.Mn 
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I. INTRODUCTION 



Recently, the discovery of superconductivity in the iron-based compounds has generated 
enormous interest, because these materials are the first non-copper superconductors with 
high superconducting (SC) critical temperature. These compounds share the same FeAs(or 
FeP) layers that are believed to be responsible for the superconductivity. Two classes of 
such compounds have been extensively investigated: (1) The LaOFeAs classes(denoted as 
FeAs-1111), space group P^/mmm, with Tc ~ 26 K through electron doping with replacing 
O- by F- fl and T c can be up to about 41 to 56 K by replacing lauthauum by other 
rare earth ions [2j, y, |J]. (2) The BaFe 2 As 2 classes(denoted as FeAs-122), space group 
I^/mmm, with Tc w 38 K through hole doping with substituting Ba 2+ for K + [5]. All 
parent_compounds show a spin-density-wave(SDW) abnormality below a temperature ~ 150 
Q] and the superconductivity is associated with the suppression of the SDW. Most 



K 



flB, 



experiment measurements have shown that these super conductors open a full gap around 
both the hole and electron Fermi pockets 0, 1(3, 
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and transport 



16, 



though the NMR Q [jjj 
171 ] experiments suggest the presence of gap nodes. At the same time, 
although the angle-resolved photoemission spectroscopy(ARPES) measurements of the band 
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221 are 



24|, a 



structure and Fermi surface(FS) in the undoped and electron doped samples 
qualitatively in agreement with the finding of the first-principle band calculation 
strong anisotropic band renormalization is found. On the other hand, the ARPES data for 
the hole doped FeAs-122 reveal that the electron FS around (n, ir) consisting of disconnecting 



patches 18|, |25j, which exhibits a significant difference from the result obtained in the 



band structure calculation. Therefore, the understanding of the role played by the electron 
correlation and the origin of the different electron FS topology between the hole-doped and 
electron-doped system is of importance. 

As the band structure calculations havin g sh own, the FS and band structures of these 



compounds are qualitatively similar 23j, |2J, |26J and the Fe-3d orbitals represent the main 



contribution to the density of states, together with a contribution from As-p orbitals, within 
several eV of the Fermi level. In the FeAs layers, the Fe atoms form a square lattice and a Fe 
atom is coordinated by four As atoms in a tetrahedron. Due to the direct Fe-Fe bonds and 
the hybridization with the As-4p orbitals, the Fe-3d orbitals form a complex band structure. 
However, the main contribution to the bands near the Fermi level comes from the d xz , d yz 



2 



and the d xy orbitals(the direct Fe-Fe bonds along the x and y axes) of the Fe atoms 



26, 



n 



281 ] . In this paper, we employ a three-orbital (the d xz , d yz and d xy orbitals) model |29|] to 



investigate the band renormalization, FS reconstruction and superconducting gap symmetry 
in the iron-based compounds with the fluctuation exchange (FLEX) approximation. We find 
that a strong anisotropic band renormalization is resulted from the Coulomb interaction 
with the strongest effect occurring around the X = (0, it) point, which is defined in the 
unfolded Brillouin zone(BZ), and this renormalization increases rapidly with the increase of 
the Hund's coupling J when J > 0.18U(U is the intra-orbital Coulomb interaction). Due 
to the band renormalization, the Fermi level for the undoped case is slightly below the 
flat band centered around the X point which is the bottom of the band along the X to 
M = (tt, 7r) direction. As a result, for the hole doped case, the Fermi level will situate below 
the flat band along the X to M direction, though it still crosses the renormalized band 
along the r to X direction. In this case, the FS around X consists of disconnecting patches. 
However, in the electron-doped case, the Fermi level is lifted to be above the flat band and 
the circular-like FS is formed. This result provides a possible explanation for the different 
FS topology observed in the electron doped and hole doped materials. We also carry out the 



same calculation based on the two-orbital model [30] , no similar FS reconstruction has been 
found. This difference is ascribed to be due to the orbital dependent renormalization. On 
the other hand, the most favored pairing state mediated by the inter-orbital spin fluctuations 
is found to be the extended s-wave with a sign change between the electron and hole Fermi 



pockets, which is consistent with the result obtained in the two-orbital model [31|,|32|,|33|. 
This indicates that the two models share the similar physics as far as the pairing symmetry 
is concerned, but exhibits difference in the band renormalization. 

The paper is organized as following. In Sec. II, we present the three-orbital model and 
introduce the FLEX method. In Sec. Ill, the numerical results for the band renormalization 
is presented and discussed. We also give a brief discussion on the pairing symmetry in this 
section. In Sec. IV, we give a summary of the results. 



II. MODEL AND FLEX METHOD 



The model Hamiltonian consists of two parts, 

H = H + H int , 
3 



(1) 



where the bare Hamiltonian Ho is given by the three-orbital model as introduced in Ref. [29]. 



In the unfolded (extended) BZ for the reduced unit cell(only one Fe atom in the unit cell as 
in Ref. 0), it can be written as H Q = ^\M k ^ k with 



Mk 



{k + Q) 

^xz,yz 

{k + Q) 

^xz,xy (k) * 
£ xz,yz{k ~\- Q) Eyg(k ~\~ Q) Eyz,xy\k) 



(2) 



y £ xz,xy(k) e *yz,xy{k) e xy{k) j 

and = (c|fj_Q, c y k z + g, c x k y ) T . Here the diagonal elements of M k denote the dispersion of 
Fe-3d orbitals d xz , d yz and d xy , while the others denote the hybridization among them. 
Keeping up to the next nearest neighbor hopping terms, we have e xz {k) = — J2 k [2ti cos k x + 
2t 2 cos k y + 4t 3 cos k x cos k y ], e yz {k) = — ^2 k [2ti cos k y + 2t 2 cos k x + 4t 3 cos k x cos k y ], e xy (k) = 



— ^ fc [2t 4 (cos k x + cos k y ) + 4t 5 cos k x cos k u 



'xz,xy 



(k) — £ *xy,xz(k) 



xz,yz (k) 
2it 7 sin fcs [29j, £yz,xy(k) 



£yz,xz(k) = -£ fc 4i 6 sinfc a sinfc„, 



'xy,yz 



(k) 



J2 k 2ii 7 sin k y [29j, 



In order to reproduce the FS and band structure feature, we set the parameters as t\ 



-1.0(« 0.4eV), t 2 = 0.7, t 3 



-0.8, t 4 



-0.3, t s = 0.2, t 6 = 0.6, * 7 



-0.35. As 



the three orbitals belong to the t 2g manifold, we set the same on-site energy to the three 
orbitals. In Fig.l, we show the band structure and FS with \i = 1.15(electron density per 
site n = 4.0) corresponding to the parent compound. We find that this three-orbital model 
can basically re pro duce the main features of the FS and band structure obtained in the LDA 



calculation 



23 



24j . The interaction between electrons is included in Hi nt as following, 



H, 



int 



2^ C Ua C lla ,C ilv' C ilv 



,<7,<t' 



+ 



i,l^l',a,a' 



+ 2 J ' Yl C lla C ila' C H'^il'a, (3) 

where U(U') is the intra-orbital (inter-orbital) Coulomb interaction, J the Hund's coupling 
and J' the inter-orbital pair hopping. 



We carry out the investigation using the FLEX approximation 34|, in which the Green's 
function and spin/charge fluctuations are determined self-consistently. For the three-orbital 
model, the Green's function G and the self-energy £ are expressed in a 3 x 3-matrix form, 




FIG. 1: (color online) Band structure and FS of the three-orbital model, (a) The band structure 
in the folded BZ with t\ = -1.0,t 2 = 0.7,i 3 = -0.8,t 4 = -0.3,t B = 0.2,t 6 = 0.6, t 7 = 0.35 and 
fj, = 1.15. (b) The FS in the unfolded BZ. The dashed lines denote the boundary of the folded BZ. 
The line with an arrow denotes the nesting vector between the hole and electron Fermi pockets. 
Panels (c) and (d) replot the hole and the electron Fermi pockets shown in Fig.(b), respectively, 
with the different colors representing the weight of the different orbitals. 
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while the susceptibility x° and the effective interaction V have a 9 x 9-matrix form. The 
Green's function satisfies the Dyson equation G(k)^ 1 = G (A;) _1 — where the self- 

energy is given by S m „(/c) = ^ J2 q V n»,mv(q)G ^(k - q) and the bare Green's function 
reads G°(k) = {\uj n — Mt + A 4 ) -1 - The fluctuation exchange interaction is given by: 

-i(f> s + f> c )x°(g)(f> s + f> c ) 

+3U S - U% m ,n„ (4) 

with spin susceptibility x s (q) — [I — and charge susceptibility x c {q) — 

[J + X (g)f/ c ] _1 X°(g)- The irreducible susceptibility is given by xj m>ni ,(g) = ~jf T,k G m(k + 
q)G mv {k). In the above, T is temperature, k = (k, \u n ) with uj n = nT(2n + 1), and / the 
identity matrix. The interaction matrix for the spin(charge) fluctuation U S (U C ) is given by: 
For i — j — k — I, U? JM = U(U^ kl = U); For i = j ? k = I, U s ijM = J{Uf jM = 2U f - J); 
For i = k,j = I and i ^ j, U s ijkl = U\Uf jkl = -U' + 2J); For i = l,j = k and i ^ j, 
Utj,ki = J{Ui jM = J); For other cases, U s ljM = 0(U^ kl = 0). 

After obtaining the renormalized Green's function G, we can solve the "Eliashberg" 
equation, 

q af3 fJLV 

*Gp v (q-k)<p^ v (k- q), (5) 
where the spin-singlet and spin-triplet pairing interactions V s and V* are given by, 

v s (q) = \u s x s (q)u s - \u c x c {q)u c + \(P S + u c ), (6) 

V\q) = -\u s X s (q)U s ~ \u c X c {q)U c + \{U S + U°)- (7) 

The most favorable SC pairing symmetry corresponds to the eigenvector (p mn (k) with the 
largest eigenvalue A. 

The Dyson equation, the self-energy and the interaction matrix Eq.(4) form a closed 
set of equations and will be solved numerically on 64 x 64 k meshes with 1024 Matsubara 
frequencies. By symmetry, we set J' = J and use the relation U = U' + 2J. In the following 
calculation, the intra-orbital Coulomb interaction U = 3.0(about 0.3 total bandwidth) is 
chosen, thus all interaction parameters are given by giving the Hund's coupling J. 

6 
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FIG. 2: (color online) Renormalization of the energy band for J = 0.2U, the red dotted(black 
solid) lines are the renormalized(bare) bands. The green, red and blue dashed lines indicate the 
Fermi levels for undoping, 10% electron doping and 20% hole doping, respectively. To show the 
renormalization more clearly, the energy bands presented here are plotted in the unfolded BZ. 

III. RESULT AND DISCUSSION 

A. Renormalization of band and Fermi surface 

We show the renormalized bands for J = 0.2U (correspondingly U' = 0.6U according to 
the relation U = U' + 2 J) together with the bare bands in Fig.2, in which the renormalized 
bands are calculated from the spectral function A(k, to) = — ^ImG^k, to) with G(k, to) the an- 
alytic continuation of the Matsubara Green's function G(k, iu) n ) by the Pade approximation. 
The Fermi levels depicted in Fig.2 are determined by calculating the number of electrons via 
the renormalized Green's function, and the undoped case is determined from the electron 
density per site n = 4.0. It is clearly seen that different bands exhibit different renormaliza- 
tion, and the strongest renormalization occurs around the X point. Compared to the bare 
Fermi level, we find that the renormailzed Fermi level is shifted up by 0.3(~ 0.12eV) for the 
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FIG. 3: (color online) (a) Renormalization factor of total bandwidth for different values of J. 
(b),(c) and (d) are the spin susceptibilities for J = 0.2C7, J = 0.05J7 and J = 0.3J7 respectively. 



undoped case. These features are consistent with the ARPES data [19|, |20[]. We define a 
total bandwidth renormalization factor as Wb/Wr, where Wb(Wr) denotes the bandwidth 
of bare(renormalized) band. The renormalization factors for different values of J are shown 
in Fig.3(a). It increases from 1.4 to 2.3 when J is increased to be around 0.5U. In particular, 
a rapid rise is observed when J is larger than 0.18U. This shows that the Hund's coupling 
plays an important role to enhance the renormalization effects. 

For U > U', the spin fluctuation will dominate over the charge fluctuation. In Fig.3(b),we 
present the static spin susceptibility x s (q, a; = 0) = EpXp, l uf( ( li a; = 0) f° r J = 0.2Z7. It 
shows four peaks around (0, ±7r) and (±7r, 0) points, which is in agreement with the neutron 
scattering experiments [fj. This arises from the nesting between the hole and the electron 
pockets connected with the vectors (0, ±7r) and (±7r, 0), as shown in Fig. 1(b). To show the 
orbital contribution to the FS, we rep lot the hole and the electron FS with different colors 
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k/n k 1% 

x X 

FIG. 4: (color online) Renormalized FS for 10% electron doping ((a) and (b)) and for 20% hole 
doping ((c) and (d)). (a) and (c) are the results shown in the unfolded BZ, while (b) and (d) in 
the folded BZ. 

representing different weights of the three orbitals in Fig. 1(c) and (d), respectively. For the 
nesting part of the FS, it is found that the d yz + d xy orbitals contribute the main weight to the 
electron pockets and the d xz orbital mainly contributes to the hole pockets. This shows that 
the spin fluctuation is mainly due to the inter-orbital particle-hole excitations. As shown 
in Eq.(3), the Hund's coupling favors the inter-orbital excitations, so the spin fluctuation 
around (0, ±7r) and (±7r, 0) will be enhanced by the increase of J. This is evidenced from 
the results shown in Fig. 2(c) and 2(d), where the spin susceptibilities for J = 0.05U and 
J = 0.3U are presented, respectively. For a small J (0.05Z7), the (0, ±7r) spin fluctuation 
disappears, while for a large J (J = 0.3U) it is enhanced. Because the band renormalization 
is increased with the rise of J as discussed above, it suggests that the strong renormalization 
around the X point mainly comes from the coupling to the spin fluctuation around (0, ±7r). 
Inspecting the details of the renormalized bands around the X point in Fig. 2, one will 
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find that the energy band near the Fermi level is renormalized strongly along the T — X 
direction, but less along the X — M direction. Another feature is that the band crossing the 
Fermi level becomes flat around the X point, and the flat portion situates at the bottom of 
the band along the V — X direction. To identify the FS after renormalization, we calculate 
the distribution of the spectral weight A(k,u) integrated over a narrow energy window 
0.02(?=s! 8meV) around the Fermi level, a method used usually in the ARPES experiments. 
The results in the unfolded BZ are shown in Fig. 4(a) for 10% electron doping and Fig.4(c) 
for 20% hole doping. To compare with experiments, their folded counterparts are shown in 
Fig. 4(b) and Fig. 4(d), respectively. In the undoped case, the Fermi level(the green dashed 
line) is just slightly below the flat band. Therefore, for the electron doped case, in which 
the Fermi level is shifted upwards as indicated by the red dashed line for the 10% electron 
doping, the Fermi level crosses the renormalized bands along both the T — X and the 
X — M directions. This gives rise to a large complete electron Fermi pocket around the X 
point as shown in Fig. 2 (a) in the unfolded BZ and in Fi g. 2( b) in the folded one, which is 



similar to th at p redicted by the LDA calculation 23|, |24|, |26( and observed in the ARPES 



experiments 19|, |21j. However, for the hole doped case, the Fermi level (the blue dashed line 
for 20% hole doping) is shifted below the flat band, so that it does not cross the energy 
band along the T — X direction anymore. Consequently, the complete electron Fermi pocket 
in the bare case is broken and the spot-like portions are formed around the X point, as 
clearly shown in Fig.4(c). After being folded, the electron FS shown in Fig. 4(d) reproduces 



the experimental results in ARPES measurements 18|, |25j]. We note that due to the band 
renormalization the crossing between the two energy band along the X — M directions is 
shifted up and approaches the Fermi level for 20% hole doping. This crossing gives rise to 
a large density of states. Therefore it adds an enhancement in the spectral weight around 
the spot-like portions and makes it be easily observed in experiments. 

The anisotropic band renormalization along the T — X and the X — M direction can be 
traced to the different orbital weights composing the energy bands along these two directions. 
Specifically, the energy band near the X point along the X — M direction is composed of 
only the d yz orbital, while that along the V — X direction is a mix of the d xy and d yz orbitals, 
as shown in Fig. 1(c) and Fig. 1(d) in which the different colors represent the weight of the 
respective orbitals. As noted above, the band renormalization is mainly due to the scattering 



off the inter-orbital spin fluctuations composed of the components Xxz,yz,Xxz,xy and Xyz 



xy 
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FIG. 5: (color online) (a) Renormalization of the energy band for the two-orbital model, the red 
dotted (black solid) lines are the renormalized(bare) bands. The green line indicate the Fermi level 
for undoping. (b) The FS for the two-orbital model, the weight of d xz and d yz orbitals are depicted 
by different colors. 



We find that the magnitude of the components Xxz,x y and Xyz,x y is the same, while it is 
larger than that of the Xxz,yz- Therefore, the d xy orbital is renormalized most strongly by 
the spin fluctuation. This suggests that the anisotropic renormalization of the energy band 
near the X point is orbital dependent. 



We note that in the two-orbital model 



30t ] the energy band near the X point along 



both directions is composed of one orbital d yz as shown in Fig. 5 (a), so the similar band 
renormalization is expected. To show this, we have carried out the same calculation for the 
two-orbital model and the result is presented in Fig.5(b). Though the strong renormalization 
still occurs in one of the bands around the X point, no strong anisotropic renormalization 
between the V — X and the X — M direction can be seen. Thus, no similar FS reconstruction 
as observed in the three-orbital model and in experiments 



18 



25] can be obtained here. 



B. Symmetry of superconducting pairing 

With the static spin susceptibility, we will further investigate the pairing symmetry me- 
diated by spin fluctuations. For J = 0.2U, which is the case with peaks around (0, ±7r) and 
(±7r, 0) points in spin fluctuations, we find that the eigenvalue A has the maximum value 
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k/n 



FIG. 6: (color online) Pairing gap functions in the unfolded BZ for U = 3.0 with T = 0.02 for 10% 
electron doping and 20% hole doping, (a), (b), (c), (d) for J = 0.2U and (e), (f) for J = 0.05U. 
(a) Ahh for electron doping, (b) A ee for electron doping, (c) A^ for hole doping, (d) A ee for hole 
doping, (e) A^h for electron doping, (f) A ee for electron doping. 



in the spin-singlet channel and it is nearly zero in the spin-triplet channel at temperature 
T = 0.02. We note that this is not the case considered in Ref. j^, where the spin-triplet 



state is expected with a larger J (such as J > U/3, see also the analysis in Ref. |33j). The 
obtained gap functions of the hole band and the electron band A ee for 10% electron 
doping and 20% hole doping are shown in Fig.3(a)-(d). It can be seen that the gap is 
basically an extended s-wave, which is nodeless around all the Fermi pockets and change 



sign between the hole pockets and t 



experiment measurements 



model 



ii 



ie electron pockets. This result is consistent with the 
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13 



181 ] and the calculation based on the two-orbital 



31] . It is noted that the pairing symmetry is similar for both the electron doped and 



the hole doped systems, although their FS topology is different as discussed above. 
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Because the spin fluctuation is dominant over the charge fluctuation, the pairing inter- 
action in the spin-singlet channel is positive (repulsive) [See Eq.(6)] and strongest around 
the wave vectors Q at which the spin fluctuation has a largest intensity. For a repulsive 
pairing symmetry, the SC gap will satisfy the condition A(k)A(k + Q) < to obtain the 
largest eigenvalue A of the "Eliashberg" equation, as can be seen from Eq.(5). For J > 0.1U, 
the inter-band (inter-orbital) spin fluctuation is dominant and has peaks around (±7r, 0) and 
(0, ±7r) which is the nesting wave vectors connecting the hole and the electron Fermi pockets. 
Thus, the gap will have an opposite sign between these two Fermi pockets. This gives rise 
to the extended s-wave. Because the SC pairing is mediated by the spin fluctuation around 
(±7r, 0) and (0, ±vr), which is the same for the three-orbital and the two-orbital models(see 
Ref. |3JJ), the same pairing symmetry will be obtained based on these two models. 




IV. CONCLUSION 

In a summary, we have investigated the band renormalization, FS reconstruction and the 
symmetry of the supercondcuting gap in iron-based superconductors using a three-orbital 
model. A strong anisotropic band renormalization due to the scattering off the inter-orbital 
spin fluctuations is found. The band renormalization is shown to be orbital dependent. As 
a result, the electron Fermi pocket exhibits different topology between the electron doped 
and the hole doped cases, which provides a natural explanation for the recent ARPES 
experiments. In addition, we have found that the most favorable superconducting pairing 
symmetry mediated by the inter-band (inter-orbital) spin fluctuations is the extended s- 
wave. 
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